Spin- and charge-density oscillations in spin chains and quantum wires 
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We analyze the spin- and charge-density oscillations near impurities in spin chains and quantum 
wires. These so-called Friedel oscillations give detailed information about the impurity and also 
about the interactions in the system. The temperature dependence of these oscillations explicitly 
shows the renormalization of backscattering and conductivity, which we analyze for a number of 
different impurity models. We are also able to analyze screening effects in one dimension. The 
relation to the Kondo effect and experimental consequences are discussed. 
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I. INTRODUCTION 

There is growing interest in impurities in low- 
dimensional electron and magnetic systems spurred by 
high temperature superconductivity and experimental 
progress in producing ever smaller electronic structures. 
There appears to be two central aspects that are stud- 
ied most in this context, namely the effect of impuri- 
ties on the transport properties in mesoscopic systems 
on the one hand,EJ and impurity-impurity interactions in 
antiferromagnetic systems due to impurity induced mag- 
netic ordercl on the other hand. In this paper we show 
that the charge- and spin-densities near impurities give 
a great deal of information about both of those aspects 
and allow us to study a number of impurity models in 
one dimension in detail. 

Induced density fluctuations at twice the Fermi wave- 
vector, so-called Friedel oscillationsjj are a common im- 
purity effect in fermionic systems, which are enhanced in 
lower dimensions. There are two distinct physical effects 
that can give rise to Friedel oscillations. The most com- 
mon source is a simple interference effect as considered 
in the original work by Friedel-u Fermions scatter off the 
impurity, resulting in a superposition of incoming and 
outgoing wave- functions. Summing up the squares of the 
corresponding wave-functions up to the sharp cutoff at 
the Fermi wave- vector fc^ results in a characteristic in- 
terference pattern with a 2kpx modulation, namely the 
Friedel oscillations. Clearly, this pattern can give a great 
deal of information about the impurity, in particular de- 
tails about the scattering process. A second source for 
the 2kpx oscillations are interaction effects due to the 
screening of an impurity with a net charge or a magnetic 
moment. A typical example of this effect is the Kondo 
screening cloud,Q which we also analyze in this paper. 
The 2kpx oscillations due to screening have typically a 
different characteristic amplitude as a function of x than 
those due to backscattering, as we will discuss in more 
detail below. 

We now consider the density oscillations in one- 



dimensional systems such as spin-chains and interact- 
ing quantum wires (Luttinger Liquids) in order to un- 
derstand the detailed effects of impurity scattering and 
screening as a function olptemperature. In the classic 
work by Kane and Fishcrtltl it was found that a generic 
impurity in a spinless Luttinger Liquid results in a renor- 
malization of the conductivity with temperature, which 
leads to a perfectly reflecting barrier at T = for repul- 
sive interactions. Interestingly, this behavior can also be 
explained in terms of repeated scattering off the Friedel 
oscillations, which gives an explicit expression oL the 
transmission coefficient in the weak coupling limit In- 
dependently, the analogous renormalization behavior was 
also found in the spin-1/2 chain,Lj where a generic per- 
turbation in the chain effectively renormalizes to an open 
boundary condition as T — ► 0. However, it is possi- 
ble that a special symmetry in the Hamiltonian reverses 
this renormalizatina, which leads to resonant tunneling 
in quantum wiresu'El og the healing of a two-link problem 
in the spin-1/2 chain.Q The renormalization behaviorJ 
that case is analogous to the two-channel Kondo effect .13' 

The renormalization flow can easily be tested numer- 
ically_bjzj examining the scaling of the finite size energy 
gapsj3't2l but we now would like to determine the reflec- 
tion coefficient directly by analyzing the induced den- 
sity oscillations which are also interesting in their own 
right. In addition, we also consider the density oscilla- 
tions from impurity models near an edge, impurities with 
a net charge or magnetic moment (Kondo-type impuri- 
ties), and integrable impurities. The detailed renormal- 
ization of the impurity backscattering as well as screening 
can be studied in each case by analyzing the induced den- 
sity oscillations as a function of temperature, which we 
determined numerically with the Transfer Ma t ri x Renor- 
malization Group (TMRG) for impurities BBeI This al- 
lows us to make predictions for conductivity measure- 
ments in quantum wires and for Knight shift measure- 
ments in spin chains, e.g. Nuclear Magnetic Resonance 
(NMR) experiments. In all cases we find a typical renor- 
malization to a fixed point of the Luttinger Liquid model, 
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which is described in terms of a simple (open or periodic) 
boundary condition in agreement with field theory calcu- 
lations. 

The rest of this paper is organized as follows. In Sec. || 
we present the model Hamiltonian and review the results 
for Friedel oscillations due to an open end (i.e. complete 
backscattering) . Different impurity models of a modified 
link, two modified links, an edge impurity, Kondo im- 
purities, and an integrable impurity are then analyzed 
in detail in Sec. III. Section IV contains a description 



of the numerical methods used and a critical discussion 
about the possible numerical errors. We conclude with a 
summary and a discussion about experimental relevance 
in Sec. [v[ 

II. THE MODEL 

The standard model we are considering here are spin- 
less interacting fermions on a one-dimensional lattice, de- 
scribed by the Hamiltonian 



H = 



E [-'(*! 



i+l 



(1) 



where n$ = is the fermion density. Although 

this Hamiltonian neglects the spin degrees of freedom 
of real electrons in quantum wires, it captures the essen- 
tial physics in conductivity experiments. Moreover, this 
model is equivalent to the spin-1/2 chain 
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where the spin operators are related to the fermion field 
by the Jordan- Wigner transformation 
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(3) 



with J = 2t, J Z = U and B = u-U. 

The model in Eq. (|l|) can be analyzed by standard 
bosonization techniques in the low-temperature limit. 
For low energies we only consider excitations around the 
Fermi-points ±kp and introduce left- and right-moving 
fermion fields with a linear dispersion relation 



- U \ _l e ikFX ip R (x). (4) 



The chiral fermion fields can then be bosonized using the 
usual bosonization rules 
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(5) 



where is the conjugate momenta to the boson field <j>. 
This results in the following boson Hamiltonian density 



(6) 



which can be solved by a simple rescaling of the bo- 
son with the interaction parameter g. The parameter 
g and the velocity v can in principle be calculated for 
any interaction strength LLand chemical potential \i with 
Bethe ansatz techniques. lij To lowest order in U we get 
.9 = 1- 2U/nv and v = y/tt 2 - n 2 + 2U/ir, so that g < 1 
for repulsive interactions. 

We now want to analyze the density oscillations using 
this formalism. Already from the decomposition of the 
fermion field in Eq. (^) it is clear that the fermion density 
may contain an oscillating component with 2kpx. To see 
this explicitly we can write the charge density in quantum 
wires (or equivalently the spin density (S z ) in spin chains) 
in terms of left- and right-movers 



RtRI 



(7) 



The first two uniform terms just represent the overall 
fermion density in the bulk system, while the last two 
"Friedel" terms are the density oscillations n osc we are in- 
terested in. In a system with translational invariance the 
left- and right-moving fields are uncorrelated (V'IV'b) = 
and no density oscillations are present. An impurity, 
however, scatters left- into right-movers and the ampli- 
tude of the oscillations gives detailed information about 
the backscattering. 

As the simplest example of this effect, let us con- 
sider an open boundary, i.e. an impurity with complete 
backscattering at the origin. In this xase-the correlation 
functions can be calculated directlyE'EJ~ril For the par- 
ticular case of the left-right correlation function at equal 
space and time we find 



{^l(x)^r(x)) oc 



ttT 



vsinh 2ttxT/v 



so that the density oscillations are given by 

irT ^ 



n osc cc sin(2kpx) 



v sinh 2nxT/v 



(8) 



(9) 



The Friedel oscillations are exponentially damped with 
temperature, because the incoming and outgoing wave- 
functions that form the interference pattern lose coher- 
ence due to temperature fluctuations. In the limit T — > 
we recover the result of Ref. ^ where a power-law decay 
of the Friedel oscillation n osc oc \/x 9 was predicted. 

It is now important to realize that the fermions or spins 
are still pinned to a lattice, i.e. x — Integer, which gives 
interesting additional effects. In particular, at half-filling 
kp = ir/2 the Friedel oscillations in Eq. ([)]) are identi- 
cally zero sin(7ra;) = for integer x, which can easily be 
understood from particle-hole symmetry (or equivalently 
spin- flip symmetry) . Half-filling is a natural state for the 
spin chains in zero magnetic field, but a small magnetic 
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field changes the Fermi vector slightly Uf = 7r/2 + B/v. 
In that case, Eq. @ becomes 



n osc oc (— l) x sin(2Bx/v) 



v sinh 2nxT/v 



(10) 



Now, the Friedel oscillations are simply alternating on 
the lattice and for distances below the magnetic length 
scale x < v/B we can use s'm(2Bx/v) — > 2Bx/v so that 
remarkably the oscillations actually increase with x 1 ^ 9 . 
This effect was first observed for the Heisenberg chain 
(J z = J, g = 1/2, v = J7r/2\y-,where the local suscepti- 
bilities xi x ) can be written asll3 

X (x)=Xo - c(-lTx hs (x), (11) 
with the amplitude of the alternating part given by 



x hs ( x ) 



T 



(12) 



V sinh AxT 

Here xo is the bulk susceptibility in the chain0 and we 
measure T in units of J. The sign was chosen so that the 
(constant) overall amplitude c of the alternating part is 
positive. The superscript bs indicates that the alternat- 
ing susceptibility is due to backscattering. As shown in 
Fig. [j] from TMRG simulations there is a characteristic 
maximum because the temperature damping eventually 
dominates over the increasing oscillations. Clearly, the 
expression in Eq. (|l2|) reproduces the shape of this al- 
ternating part rather well, although we have neglected 
possible logarithmic corrections (multiplicative and ad- 
ditive), which may be responsible for the apparent shift 
in the characteristic maximum in Fig. [l]. The numer- 
ical TMRG results of the local susceptibility near the 
open end x bs ( x ) wm be used as the reference data for 
a completely backscattering impurity in our studies in 
the next section. The numerical data automatically con- 
tains all corrections due to irrelevant higher order op- 
erators. The logarithmic corrections to Eq. (|lj) due to 
the leading irrelevant operator have a special behavior 
near a boundary,EU which we have not tried to predict 
for the local susceptibility, but numerically we find that 
a possible multiplicative logarithmic correction for x( x ) 
appears to have a negative power of ln(a;). The max- 
imum in Fig. [l] occurs at x oc 1/T with an amplitude 
Xalt °c which results in a characteristic feature 

in NMR experiments, so that it was possible to confirm 
this effect experimentally as well.ES At zero temperature 
the ground state has a staggered magnetization which 
has a maximum in the copter of a finite chain (assuming 
an odd number of sites). c3 The magnetization for finite 
chains with impurities has also recently been analyzed, 
which resulted in interesting patterns that repeal the na- 
ture of the strong correlations in the systemO 

Even for a partially reflecting impurity we expect that 
the same alternating contribution as in Eq. (|12| ) due to 
backscattering is present, but with an amplitude c that 



increases monotonically with the reflection coefficient R. 
In fact we can make a firm connection between the rel- 
ative amplitudes and the reflection coefficients by con- 
sidering free fermions U = for which we can find the 
eigenf unctions exactly even in the presence of impurities. 
Clearly the eigenfunctions are given by plane wave solu- 
tions |fc) which contain a special mix of left- and right- 
moving components due to the impurity. Just like with- 
out impurities there are in fact always two such degen- 
erate orthogonal solutions. We found the solutions for 
generic impurity models and looked at the spatial struc- 
ture of the square of the wave- functions, which contains 
an interference pattern of incoming and outgoing waves. 
In general, we always find 
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|(a;|fc)| 2 = - ( 1 + -vAR(fc) cos(2fcx + 2$) 

7T V 



(13) 



where the summation over the two degenerate solutions 
is implied. Here R(k) is the ordinary k-dependent reflec- 
tion coefficient which has been determined independently 
according to text-book methods. Therefore, the magni- 
tude of the interference is exactly given by the square 
root of the reflection coefficient, which is maybe not too 
surprising but very useful in our analysis. In particular, 
when we consider the fermion density at half filling we 
are really directly looking at the spatial structure of the 
wave- function. We can write near half-filling (i.e. for a 
small field B in the spin chain model) 

n(x)-l/2 = / \(x\k)\ 2 dk 

J it/ 2 



fl =°f \m\ 2 

= B [xo - c R {-lf X hs ] 



(14) 



where we have used the fact that the spin density for 
the Heisenberg chain in a small field is just given by the 
susceptibility in Eq. ([ll]) , but with a coefficient cr which 
now depends on the reflection coefficient R near half fill- 
ing. Together with Eq. (13) we therefore arrive at the 
central result that at half-filling the reflection coefficient 
is proportional to to the square of the alternating density 
amplitude 



R 



(15) 



where c = Cr=i is the coefficient corresponding to com- 
plete backscattering in Eq. (pd|). We use this formula to 
estimate the reflection coefficient from the density oscil- 
lations for various impurity models in the following. 
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FIG. 1. Local susceptibility close to the open end of a 
spin- 1/2 chain from TMRG data for T = 0.04 J compared 
to Eqs. (0) and (|l|) with c = 0.51, which was determined 
by matching the characteristic maxima. 

As mentioned above there may also be 2kpx density 
oscillations due to screening, so that the alternating sus- 
ceptibility is in general a sum of two parts 

X aU (x) ee X (x) - X o= (~IT [ X scrooning (^) - c RX hs (x)} . 

(16) 

In the case of overscreening the neighboring spins (or 
electrons) overcompensate the magnetic (or electric) im- 
purity and leave an effective impurity with opposite mo- 
ment which in turn gets screened by the next nearest 
neighbors and so on. This finally results in a screening 
cloud. Screening is purely an interaction effect where a 
2k px density oscillation is induced by an "active" impu- 
rity Hamiltonian (x/j' L x/j R Hi mp ) ^ 0. The 2kpx oscilla- 
tions due to backscattering, however, are purely an in- 
terference effect and are even present in non-interacting 
fermion systems. The special shape and the increasing 
nature of the alternating part in Eq. ( p"2| ) for g = 1/2 
makes it possible to easily identify the contribution due 
to backscattering, so that we can always separate the two 
possible effects near half-filling. In what follows we there- 
fore always use the special choice of coupling U = 2t cor- 
responding to the Heisenberg model J z = J. This model 
can be used to demonstrate the generic behavior of im- 
purity effects in mesoscopic systems and also gives ex- 
perimental consequences for spin-chain compounds. The 
Luttinger Liquid parameter takes the value g — 1/2 in 
this case, which is the strongest possible interaction at 
half- filling before Umklapp scattering becomes relevant. 



III. IMPURITY MODELS 

A. One modified link 

Maybe the simplest impurity to consider is a weak link 
in the chain, i.e. a modified hopping J' between two sites 
in the chain as shown in Fig. |2| 



i+1 



t+1 % 



(17) 
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FIG. 2. One modified link. 

The wave- functions and reflection coefficient R(k) for 
this problem can be calculated exactly, with the result 
that 



R(k) = 



r - 2t z j' 



2 t/2 



J" 



t 4 - 2t 2 J' 2 cos 2k + J' 4 ' 



(18) 



However, once the interaction U is introduced this prob- 
lem becomes highly non-trivial and the reflection coeffi- 
cient renormalizes with temperature T . The interacting 
system has been studied in the context of both spinless 
fermionsEI and the spin-1/2 chain,!] where it was found 
that repulsive interactions U > make the perturbation 
of one link relevant, so that it renormalizes to a com- 
pletely reflecting barrier as T — > 0. A small weakening of 
a link J' < t produces a relevant backscattering operator 
in the periodic chain of scaling dimension d = <?, so that 
this link effectively weakens further as the temperature 
is lowered. Below a cross-over temperature Tk (analo- 
gous to a Kondo-temperature) the link has weakened so 
much that it is more useful to consider the problem of 
two open ends that are weakly coupled, which is now 
described by an irrelevant operator of scaling dimension 
d = l/g. Therefore, this coupling weakens further and 
ultimately the open boundary condition represents the 
stable fixed point as T — > 0. The same analysis is also 
true for a slight strengthening of a link J' > t, because 
in this case the two ends lock into a "singlet" state as 
the effective coupling grows, and the remaining ends are 
weakly coupled with a virtual coupling of order t 2 / J' 
which is again irrelevant. 

We consider the interacting system with U — 2t, which 
we can write in terms of an SU(2) invariant spin Hamil- 
tonian via the Jordan- Wigner transformation in Eq. (|^) 
with a modified Heisenberg coupling between two spins 



H = J^S, 

i#0 



J So • Si. 



(19) 



We now want to analyze the density oscillation near the 
impurity in order to extract the reflection coefficient as 
described above. In Fig. || we show the amplitude of the 
alternating spin density for different coupling strengths 
J'. Clearly the shape as a function of distance x remains 
largely the same as in Fig.[l| for all J' so that the func- 
tional dependence in Eq. ( |X2[ ) is still adequate, but with 
an overall coefficient c which is now related to the reflec- 
tion coefficient R as postulated in Eq. (|i~5|). 
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FIG. 3. Envelope of the alternating susceptibility of the 
one- link impurity at T = 0.04J for J'/ J = 0.0, 0.2, 0.4, 0.6, 0.8 
from above. 

The reflection coefficient is directly related to the 
renormalization behavior above. The basic idea behind 
renormalization is to use an effective Hamiltonian with 
renormalized parameters as a function of T . To esti- 
mate the reflection coefficient it is therefore possible to 
make a simplified but intuitive analysis by using the free 
fermion result in Eq. (|l8|), but with a renormalized cou- 
pling strength J'(T). Below the cross-over temperature 
T < Tk, the effective potential is small and given by the 
renormalization behavior of the leading irrelevant opera- 
tor J'(T) cx J'T 1 ^- 1 . This results in 

1 - Roz j ,2 T 2/g - 2 , (20) 

which is the universal behavior near the stable fixed point 
as first predicted in Ref. [j]. Above the cross-over tem- 
perature T > Tk the renormalization behavior is better 
described by a relevant operator on the periodic chain 
giving J - j'(T) cx (J - J^T 3 ' 1 . From this result it 
would even seem that we can recover the periodic chain in 
the high temperature limit, but it is of course important 
to realize that the renormalization is no longer possible 
above a cutoff of order J. For an initial bare coupling 
J' ~ J very close to the unstable fixed point Tk <C J we 
therefore find that the effective coupling stops renormal- 
izing at its bare value J' — > J' for large T. In summary, 
the temperature dependence above Tk is not as universal 
as m Eq. @, but we may still write 

R(x{J-J') 2 , (21) 

for J' - J and T > T K . 

It is now straightforward to extract the relative coeffi- 
cient cr/c in Eq. ( |i~4| ) from the numerical data by simply 
dividing the amplitude of the alternating part for each 
coupling J' in Fig. || by the reference data of x bs for the 
open chain. According to Eq. ( |l5|) the square of this rela- 
tive coefficient then gives the reflection coefficient. Fig. |] 
shows the results for the temperature dependent reflec- 
tion coefficient from our TMRG data. The renormaliza- 
tion to a perfectly reflective barrier can clearly be seen 



as T — > 0. The behavior for couplings close to the peri- 
odic fixed point (J' > 0.4J) is consistent with Eq. (|2l|). 
For smaller couplings the cross-over temperature Tk is 
larger, and we see an extended region where the scaling 
of the stable fixed point with J' 2 and T 2 /^ 2 in Eq. @) 
holds (here g = 1/2). We can also compare our results 
to the findings of Matveev et al in Ref. ^ where an ex- 
plicit formula for the transmission coefficient was given 
l-Rot [(D/T) 2a R /(l-R Q ) + I]- 1 in terms of the non- 
interacting reflection coefficient Rg in Eq. (|i~8|), a cut-off 
D, and a small interaction parameter a = 1/g — 1. Un- 
fortunately, the interaction parameter is large in our case 
a = 1 so that this formula does not quantitatively agree 
with our findings in Fig. ^. Qualitatively, their result 
looks rather similar, but we observe a sharper renormal- 
ization at low temperatures near the unstable fixed point 
(J' > 0.4 J). Indeed we find that the region where the 
famous scaling in Eq. (|2(]) is valid turns out to be ex- 
tremely narrow for J' > 0.4J. 

Another aspect is the high temperature behavior where 
the non-interacting reflection coefficient in Eq. JTs| ) 
should be approachedB. This is indeed the case near the 
unstable fixed point J' > 0.4J where the non-interacting 
value is quickly reached with high accuracy. However, 
near the stable fixed point (J' < 0.4J) we find that the 
reflection coefficient can renormalize even well below the 
non-interacting value, so that the interactions actually 
enhance the conductivity at higher temperatures in this 
case. The reason for this unexpected behavior is that 
the cross-over temperature is larger than the cut-off near 
the stable fixed point Tk 3> J, so that the renormaliza- 
tion may continue beyond the bare coupling constants at 
higher temperatures. 
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FIG. 4. Reflection coefficient R of one modified link for 
J' I J = 0.1, 0.2, 0.4, 0.6, 0.8 from above. The lines are guides 
for the eye. 

B. Two modified links 

We now consider the impurity of two neighboring mod- 
ified links in the chain as shown in Fig. 0. For the in- 
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teracting case U — 2t we can again write this model in 
terms of a Heisenberg spin chain model 



H = J 



E 

i#-1.0 



Si • S,+i + J'Sq ■ (S_i + Si) 



(22) 



This type of impurity may correspond to a charge island 
that is weakly coupled to a mesoscopic wire or to doping 
in a quasi-one dimensional compound where one atom in 
the chain has been substituted. We have recently con- 
sidered this type of impurity in the context of doping in 
spin-1/2 compounds and as a simple .experimental exam- 
ple of the two channel Kondo effect B In this section we 
analyze the induced density oscillations in more detail, 
especially in connection with the reflection coefficient. 
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FIG. 5. Two modified links. 



The model in Eq. (|22j) is equally simple as the one- 
link impurity, but thegenormalization behavior is known 
to be quite different. Already for the non- interacting 
case at half filling the system shows a resonant behavior 
with perfect transmission R = 0, so that this corresponds 
to the simplest case of resonant tunneling considered by 
Kane and Fishema (at half-filling the impurity potential 
is automatically tuned to the resonant condition). With 
interactions U ^ the reflection coefficient is no longer 
exactly zero, but shows nontheless a renormalization to 
perfect transmission as T — > in sharp contrast to the 
one-link impurity. This difference in renormalization be- 
havior is easily explained by the different parity symme- 
try of the problem (namely site- instead of link-parity) . 
For a small perturbation from a periodic chain J' ~ J 
the leading operator is now irrelevant with scaling dimen- 
sion of d = 1 + <?, so that a perfectly transmitting chain is 
the stable fixed point. For small couplings J' > on the 
other hand, the leading perturbing operator is marginally 
relevant, and the situation is similar to the two channel 
Kondo effect where the two ends-,ftf the chain play the 
role of two independent channels.Lrcl 

Apart from the renormalization behavior there is an- 
other key difference between the one- and two-link impu- 
rities: In the two-link impurity model there is an "active" 
impurity site that carries a spin or charge degree of free- 
dom, which in turn must be screened by the surrounding 
system. Therefore, the density oscillations are no longer 
simply determined by the backscattering in Eq. ([l2]) , but 
there is also a so-called screening cloud induced in the 
system. From perturbation theory in the leading irrele- 
vant operator the functional dependence of this screening 
cloud can be calculated^ and the total alternating density 



!lt is a sum of two contributions 



X al \x) = di-lf ln[coth(arT)] - c H (-l) 3 'x bs (x), (23) 

where the first term is the induced screening cloud while 
the second term is the familiar contribution due to 



backscattering in Eq. (|l2|). Interestingly, the two con- 
tributions have opposite sign, so that the density oscilla- 
tions vanish at a special distance from the impurity, but 
then increase again due to the backscattering contribu- 
tion. This behavior is shown in F ig. H together with a 
fit to the two contributions in Eq. (J230T The special dis- 
tance at which the density oscillations vanish grows as 
we approach the stable fixed point (J' — > J or T — > 0). 
As already with the one-link problem, we use again the 
numerical open chain data as a reference for x hs instead 
of the more simplified analytical form of the backscat- 
tering contribution in Eq. (12) since this minimizes the 
corrections due to irrelevant operators. However, even 
the analytical form in Eq. ( p^ ) gives very good fits, so 
that none of our our findings are affected by this choice. 
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FIG. 6. Alternating part of the local susceptibility for the 
two-link impurity for T/J = 0.04 and J'/J = 0.6. Fit to 
Eq. (§|. 

It is now straightforward to extract the reflection coef- 
ficient from the numerical data with the help of Eq. ( |l5| ) 
and Eq. ( f23| ) as shown in Fig. [?]. Below a cross-over 
temperature Tk depending on J' the reflection coeffi- 
cient clearly decreases and eventually approaches perfect 
transmission as T — > 0. Above Tjc the renormalization of 
the reflection coefficient is rather weak and converges to 
a finite constant (never approaching complete reflection 
as the temperature increases). 
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FIG. 7. Reflection coefficient R of the two-link impurity 
for J' I J = 0.05,0.1,0.2,0.4,0.6,0.8 from above. The lines 
are guides for the eye. 

Equally interesting is the induced screening cloud. In 
this case, the coefficient cj approaches a constant as 
T < Tk as it should, since this contribution was deter- 
mined from perturbation theory around the stable fixed 
point. Above the cross-over temperature, however, this 
contribution vanishes quickly. This behavior is shown 
in Fig. ^: In general the behavior of the coefficient cj 
vs. J' is temperature dependent and ci increases as the 
temperature is lowered. However, as T <C Tk all curves 
approach a limiting value, which gives a universal behav- 
ior as a function of J' (thick line). 



havior of the two channel Kondo fixed point. Note, that 
although the coefficient cj is finite as T — > 0, the screen- 
ing cloud itself diverges logarithmically with — ln(xT), 
which is a clear indication of the famous over-screening 
in the two channel Kondo effect. As we approach the 
unstable fixed point the order of limits becomes cru- 
cial: For zero coupling there is no screening cloud at 
all limT_>o limj'^o c i = 0, while for zero temperature 
the coefficient becomes infinite limj/_»o lim-r^o cj = oo. 
Remarkably, exactly at zero temperature a minute per- 
turbation therefore induces an infinite screening cloud, 
although this behavior occurs in an unphysical limit. 

C. Impurity at the edge 

Another category of impurities we can consider are im- 
perfections near the end of a chain. In this case the 
boundary always gives complete backscattering, but as 
we will see the impurity can still give interesting effects 
on the density oscillations. The simplest case to consider 
is a modified link at the edge of a chain as depicted in 
Fig. |j| For the interacting case U = It it is again use- 
ful to write the Hamiltonian in terms of the Heisenberg 
spin-chain model 

oo 

H= J^S i -S i+1 + J'So-Si. (24) 
»=i 
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FIG. 8. Coefficient cj vs J' of the two link impurity for 
different temperatures T/J = 0.2, 0.1, 0.04, 0.025, 0.0167, 0.01 
from below. For J' ~ J and/or low temperatures ci ap- 
proaches a universal T- independent curve (thick line). Inset: 
a vs T. The lines are guides for the eye. 

The competing contributions in Eq. ( p3|) have the op- 
posite renormalization behavior: Above Tk backscat- 
tering is constant, while the screening cloud is reduced 
which is the open chain behavior. Below Tk on the other 
hand backscattering is reduced, while the coefficient for 
the induced screening cloud is constant, which is the be- 



J' 
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FIG. 9. Edge impurity. 

Just like the two- link impurity was related to the two- 
channel Kondo problem, we can identify the field theory 
description of the edge impurity model with the regu- 
lar one-channel Kondo problem. There are two possible 
fixed points: The case J' = corresponds to the unsta- 
ble fixed point of a decoupled spin at the end of a chain 
with a marginally relevant perturbation for J' 0. The 
case J' = J corresponds to the completely screened spin, 
which is a stable fixed point with a leading irrelevant 
operator of scaling dimension d = 2. Just like in the or- 
dinary Kondo effect both fixed points are represented by 
the same boundary condition and differ only by a simple 
7r/2 phase shift on the fermions. (The infinite coupling 
fixed point J' — > oo is also stable, but is actually abso- 
lutely equivalent to the J' = J fixed point since both 
cases represent a ir/2 phase shift on the fermions by re- 
moving or adding a site, respectively). For intermediate 
couplings the phase shift <& takes on values between and 
7r/2 which will be reflected in the backscattering contri- 
bution of the density oscillations as we will see below. 

A screening cloud for the impurity spin at the end 
should also be present in this model, but with a different 
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behavior than for the overscreened case in Eq. (p3[). In- 
stead we find that the leading operator that causes the 
screening cloud is the same as that for an edge magnetic 
field in the xxz-chain which has been analyzed in Ref. |2^, 
so we can use the corresponding result for the shape of 
the induced screening cloud. Taking into account finite 
temperatures and the phase shift on the fermions we can 
write for the density oscillations 



X alt (x)=cj 



{-i) x VT 

^sinh(4xT) 



cos(ttx + 2$) c x (x), (25) 



where the first term is the induced screening cloud, while 
the second term is the backscattering contribution in 
Eq. (12) but with a phase shift $. However, the coeffi- 
cient c always takes the value corresponding to complete 
backscattering in Eq. ([ll]). There is also an implied shift 
of 2$ /V in the argument of x bs , which we used for a self- 
consistent fitting. The effective boundary condition in 
the continuum limit is therefore technically between two 
lattice sites (although it is not really that meaningful to 
define locations on the scale of less than a lattice spacing 
in the continuum limit theory anyway). 

Figure |l^ shows the envelope of the alternating part of 
the susceptibility for temperature T = 0.04 J and differ- 
ent couplings J', which always fits well to the superposi- 
tion in Eq. (p5|). At the fixed points J' — and J' = J 
there is no screening, but the backscattering contribution 
has opposite signs due to the 7r/2 phase shift. 




FIG. 10. Alternating susceptibility for the edge impurity 
at T/J = 0.04 for J' /J = 0,0.1,0.2,0.4,0.6, 1.0 from above. 
Fits to Eq. (||). 

It is now straightforward to extract the screening cloud 
amplitude c/ and the phase shift <f> from our numerical 
data for all temperatures and couplings J'. As expected 
we find that the phase shift increases with J' and renor- 
malizes to lar ger values as the temperature is lowered as 
shown in Fig. O. In the limit of very low temperatures 
the jump to the stable fixed point value <E> = ir/2 becomes 
more abrupt as a function of J'. 



<F 




FIG. 11. Phase shift of alternating part of the edge impu- 
rity for T/J = 0.25, 0.1, 0.04, 0.02, 0.01 from below. The lines 
are guides for the eye. 

The screening cloud coefficient ci again approaches a 
constant as we lower the temperature below Tk as shown 
in Fig. |l2|. Although formally the behavior looks similar 
to the over-screened case of the two link problem in Fig. ^ 
it is important to realize that now the screening cloud in 
Eq. ( p5| ) is finite as T — ► and drops off with 1 jx (while 
in the two link case the screening cloud was divergent 
with hixT). 
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FIG. 12. Coefficient ci vs. J' of the edge impurity for 
T/J = 0.2,0.1,0.04,0.02,0.133,0.01 from below. The lines 
are a guide for the eye. 



D. Generalized two link impurity 

It is now instructive to summarize the findings of the 
three impurity models in the previous subsections by con- 
sidering one generalized two link impurity model that is 
not symmetric as shown in Fig. [f| 

H = J J2 s i- s i+i + AS_i-S + J 2 S -Si. (26) 

i 5^-1,0 



The three impurity cases above can be identified easily: 
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•h ^ Ji — J one modified link in Eq. (19) 



• Ji = J-2 7^ J two modified links in Eq. ( |22|) 

• Ji = 0, J 2 ^ J edge impurity in Eq. ( |24| ) 

The density oscillations for the more general model in 
Eq. ( p6| ) are much more complex than in the special cases, 
so that a detailed analysis of this effect is not always use- 
ful. The renormalization behavior on the other hand is 
straightforward and can be read off from what we already 
know about the special cases. 

J, L 



• • • » • • • — 

-3-2 -10 12 3 
FIG. 13. Generalized two-link impurity. 

A weak coupling J\ > and J2 > to an additional 
site is always marginally relevant, so that the open chain 
with a decoupled impurity site is unstable for any antifcr- 
romagnetic coupling (i.e. negative hopping probability). 
The periodic chain on the other hand is only stable for 
the special site-parity symmetric case J\ = J2 , where the 
renormalization behavior is analogous to the two channel 
Kondo effect. In general, however, one of the two cou- 
plings is larger and renormalizcs to unity, absorbing the 
spin. The smaller coupling is then irrelevant as in the 
one-weak problem, so that the stable fixed point is an 
open chain with an absorbed impurity site J\ = J, J2 = 
(or J2 — J, J\ — 0) in most cases, except for a site-parity 
symmetric impurity or two ferromagnetic coupling con- 
stants. The complete renormalization flow is summarized 
in Fig. [IJ where the possible fixed points are indicated by 
the black dots. In cases where the coupling diverges to in- 
finity a singlet forms, and we can therefore again describe 
the system by one of the four finite fixed points in the 
figure. Interestingly, the more stable fixed points always 
have a lower ground state degeneracy, in accordance with 
the g-theorem.cj The phase diagram in Fig. [lj] is valid 
for all interaction strengths < U < 2t as long as the 
system is half-filled. 



\ ' 
1 




\/0' 


> < 

y 



t « HI 


' J 




v- *- « 

, ^ 


> — ■* > 

' v S 



FIG. 14. Renormalization flow diagram. 



E. Spin-1 impurity 

We now turn to a magnetic impurity in the chain with 
spin 5i mp = 1 given by the Hciscnbcrg Hamiltonian 



H 



Si+l + J'Si, np • (Sq + Si) 



(27) 



as shown in Fig. |l5|. In the previous impurity models 
in Sections [II A - III E)| it was always possible to inter- 
pret the Heisenberg Hamiltonians equally well in terms 
of mesoscopic systems and electrons hopping on the lat- 
tice by identifying the spin-1/2 impurity in terms of an 
extra site or charge island. However, for the spin-1 im- 
purity in Eq. ( p7|) no meaningful interpretation in terms 
of spinless fermions is possible. On the other hand this 
impurity model has important implications for doping in 
quasi one-dimensional spin-1/2 compounds, so that we 
find it useful to discuss it here. 



J'/"\ J' 
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FIG. 15. The spin-1 impurity. 



Similar to the impurity models in Sections III B and 



[II C we find again that the field theory language is anal- 
ogous to a Kondo impurity model. The two ends of the 
spin-chain play the role of the two channels coupled to 
a spin-1 impurity. A small antiferromagnetic coupling 
is therefore marginally relevant and the renormalization 
flow goes to the strong coupling limit. The stable fixed 
point is given by an open spin chain with two sites re- 
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moved and a decoupled singlet containing the spin-1 and 
the two end spins {J 1 — > oo). 



20 



Just like the edge impurity in Sec. IIIC this Kondo- 
type model is an exactly screened impurity. The shape 
of the screening, cloud is again given by that of an edge 
magnetic fieldcl just like in Eq. (|25| ) 



X M (x) = cj 



v /sinh(4a;T) 



cb(-1)V(z), 



(28) 



where the first term is again the induced screening cloud, 
while the second term is the backscattering contribution 
in Eq. (|l2]). As shown in Fig. [l^ the fits to this expres- 
sion are excellent (again using the open chain data as 
a reference for x bs )- The coefficient cj for the induced 
screening cloud again approaches a constant for temper- 
atures below Tk which results in a universal curve as 
T — > as shown in Fig. The backscattering coef- 
ficient is an indication of the effective phase shift and 
changes sign depending on the temperature and coupling 
strength. From Fig. [16| it is clear that the backscatter- 
ing coefficient cr is positive for small coupling strengths 
J' (or equivalently high temperatures) and negative for 
larger coupling strengths J' (or equivalently lower tem- 
peratures). The renormalization of cr is explicitly shown 
in the inset of Fig. |TJ. As T — ► the jump of cr to 
negative values happens at smaller J' and becomes very 
sharp. 




FIG. 16. Envelope of alternating part at T/J 
for the spin-1 impurity. From above at x 
J'/J = 0, 0.05, 0.1, 0.2, 0.4, 0.8, 4.0. Fits to Eq. 
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FIG. 17. Coefficient a of the spin-1 impurity for T/J — 
0.2,0.1,0.04,0.025,0.0167,0.0133,0.01 from below. Inset: 
Backscattering coefficient cr. The dashed lines are a guide 
for the eye. 

More interesting are the experimental consequences for 
Knight shift experiments in doped spin-1/2 chain com- 
pounds (as for example Ni doping in CuO chains). For 
that case we can predict an interesting NMR spectrum 
with a characteristic feature (sharp edge) corresponding 
to the maximum in the alternating susceptibility. Such 
a sharp edge has been observed before in NMR experi- 
ments oh. spin-1/2 chain compounds with non-magnetic 
defects B3 In that case the sharp edge broadens with a 
1/VT behavior as discussed in Sec. |l[ For the mag- 
netic spin-1 impurities a sharp edge from the maximum 
in the backscattering part may also be present, but it 
depends on if the temperature is above or below Tk how 
this feature changes. Above Tk the backscattering part 
becomes weaker as the temperature is lowered, but the 
induced screening cloud increases, so that the sharp kink 
may vanish in a quickly broadening line-shape from the 
screening cloud as shown in the left part of Fig. |l8|. Be- 
low Tk on the other hand, the screening has saturated 
and the backscattering contribution dominates again (al- 
beit with a phase shift). Therefore, the kink feature in 
the NMR spectrum will sharpen further as the tempera- 
ture is lowered and widen with the usual 1/yT behavior 
as shown in the right part of Fig. [l|. The detailed T- 
dependence can be predicted for any particular value of 
J' of an actual experimental compound. 
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FIG. 18. NMR signal of the spin-1 impurity for J'/ J = 0.1 
(left) and J'/J = 1.4 (right). 



F. Integrable impurity model 

Finally, we would like to consider a more exotic im- 
purity model which has been especially constructed to 
preserve the integrability of the entire system.EH We con- 
sider here the simplest non-trivial example of such an 
impurity model which corresponds to an impurity spin 
with Si mp = 1 that is coupled in a special way to two 
sites in the chain. The corresponding Hamiltonian has 
been set up in Ref. |2?] 



H — J Sj ■ Si + i — "TfSo ' Si 

+ 4f [(So + Si) • S; mp + {So • Si mp , Si • Simp}] > 



(29) 



where Si mp is the external spin-1 impurity and {,} de- 
notes the anticommutator. 
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FIG. 19. The integrable impurity. 

A closer analysis of this modelil showed that the ther- 
modynamics at low temperatures were in fact described 
by a periodic spin chain with one additional site and an 
asymptotically free impurity spin with S = 1/2, so that 
it appears that the original spin-1 has somehow been par- 
tially absorbed by the chain. From a field theory point 
of view it was later shown that this type of impurity cor- 
responds in fact to an unstable fixed point which can 
only be reached by an artificial tuning of the coupling 
parametersEj 

We are now interested in what kind of density oscil- 
lations might be observable from such an impurity. In- 
terestingly, we found that the density oscillations were 
identically zero at all temperatures as if the system was 
translationally invariant. The impurity Hamiltonian in 



Eq. (g^) was of course constructed in a way to avoid all 
backscattering, but it is remarkable that even the induced 
alternating part from the magnetic impurity vanishes ex- 
actly, i.e. no conventional screening takes place. 

Nonetheless, the impurity spin is somehow reduced 
from a spin-1 to an effective spin- 1/2 as the temperature 
is lowered. This can be explicitly seen from the impurity 
susceptibility in small magnetic fields 



{Si 



imp/ 



B 



a 



Curie 



T 



(30) 



where we have assumed some type of Curie-law. At high 
temperatures the impurity susceptibility must follow the 
Curie-law for a spin-1 Ccuric = 2/3, while at low temper- 
atures a Curie-law for a spin- 1/2 Ccuricrr 1/4 has been 
predicted up to logarithmic corrections SB In Fig. ^ we 
plot the temperature dependent Curie constant (i.e. the 
impurity susceptibility times temperature). It appears 
that the asymptotic value Ccuric = 1/4 is indeed ap- 
proached with logarithmic corrections as T — > 0. The fit 
in the figure is 



a 



Curie — T 



1 



81n(27r/T) 
with a = 1.62 and b = 1.32. 



0.45 



In Qn(2n/T)/b) 
1 ln(27r/T) 2 



(31) 
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FIG. 20. Susceptibility of the spin-1 in the integrable model 
multiplied by temperature. Fit to Eq. (f$l|). 



IV. NUMERICAL METHOD 

The numerical method we have used here is based an 
the Density Matrix Renormalization Group (DMRG)cJ 
applied to transfer matrices. While the ordinary DMRG 
considers the properties of individual eigenstates in a 
finite system, we are interested in the thermodynamic 
limit, namely properties of an infinite system at finite 
temperatures. This can be achieved by -the Transfer 
Matrix Renormalization Group X,T-MRG)J13 which we 
adapted especially for impuritiesoila as we will review 
briefly. We consider the partition function Z of the mod- 
els in Eqs. (|l|) and (|2|). After the standard Trotter de- 
composition, we obtain for an infinite system [L — > oo) 
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Z = lim trT 



L/2 



M- 



,1/ 



lim A 

Af- 



L/2 
M i 



(32) 



(S?)os C ex (1 - e )* = exp(-je) 



(35) 



where Tm is the transfer matrix with M time-slices. In 
the limit of infinite system size only the largest eigenvalue 
Am determines the thermodynamics of the system, which 
we find numerically. We start with small time-steps so 
that the Trotter-error is negligible, and successively in- 
crease the number of time-slices M to reach lower tem- 
peratures. At each step the dimension of Tm increases 
so we keep only the most important states to describe 
the state with the highest eigenvalue Am by using the 
DMRG algorithm with some modifications for asymmet- 
ric matricesJlJ A measurement of the local spin-density 
at site j for example is straightforward, since we can just 
absorb the measuring operator Sj into one of the transfer 



matrices T-. 



M 



T 



M 



1 

z 



(5?) = 4 tr 5? e~ m 



A 



M 



(33) 



where (4>m\ and \i/>m) are the left and right target states 
for the eigenvalue Am- So far we have considered a trans- 
lational invariant system. 

We now introduce a generic impurity which modifies 
one of the transfer matrices Tm — > T; mp . Even in the 
presence of impurities the thermodynamics of the sys- 
tem is entirely determined by the highest eigenvalue Am 
and corresponding eigenstate of the pure transfer ma- 
trix Tm which always appears with an infinite power in 
the partition function in Eq. d32|). The measurement of 
the spin (or charge) density near the impurity is again 
straightforward. For the spin density at a distance of j 
sites from the impurity we write 



(il> M \Tft {T M y/ 2 T imp |^M) 



A 



i/2+1 

M 



(lpM\T imp \4>M} 



(34) 



Since we step-wise approximate the transfer matrix, 
it is important to make a careful error-analysis. The 
error due to the Trotter approximation is the simplest 
to estimate since it is just proportional to the square 
of the time-step r = l/TM. We found that a value of 
t = 0.05/J makes this error negligible compared to the 
DMRG truncation errors. To estimate the truncation 
errors we can compare our results to the exact solution 
of the free fermion Hamiltonian in Eq. (Q) with U = 0. 
The structure of the transfer matrix is not fundamentally 
changed by taking U — so that the truncation error will 
be of the same order as for U ^ 0. Keeping 64 states we 
find for the local response of the spins closest to typical 
impurities a relative error of less than 10 -4 for T > 0.04, 
less than 10~ 3 for 0.02 < T < 0.04 and a relative error of 
less than 10 -2 for temperatures 0.01 < T < 0.02. How- 
ever, already from Eq. ( |34| ) it is clear that the spin and 
charge densities far away from the impurity will contain a 
larger error. Each transfer matrix contains a small error 
e which then gets exponentiated in Eq. ([54j) and hence 
the oscillating part of the density (S*) is suppressed ex- 
ponentially with distance j 



where e depends only on temperature. This exponential 
suppression with the distance from the boundary is again 
a consequence of the fact that the incoming and outgo- 
ing waves lose coherence but this time due to error fluc- 
tuations. However, the corresponding energy scale from 
the truncation error is always smaller than the tempera- 
ture in our case. We observe that the suppression error 
in Eq. (35) is actually very systematic, so that we can 
even correct our data very well using Eq. (|35|). For free 
fermions we find to high accuracy the following depen- 
dence of the error on temperature 



e = 0.06cxp(-58T), 



(36) 



where we have kept 64 states in the TMRG simula- 
tions. For interacting fermions the suppression also has 
the exponential dependence in Eq. (|35|), but the en- 
ergy scale e is in general dependent on the interaction 
U. For the Heisenberg model an independent analysis 
of the free energy hinted at a value of approximately 
e = 0.02 exp(— 34T), but the value in Eq. (36) is more 



reliable and gives a relatively good estimate of the error 
for all interaction strengths. We chose to correct our data 
for the alternating fermion densities by dividing out the 
factor in Eq. ( p5| ) together with the estimate in Eq. ( |3"6| ) 
in all cases presented above. However, the use of this 
correction or the particular choice of the error e makes 
no qualitative difference in any of our findings, since the 
temperature suppression always dominates (i.e. the en- 
ergy scale in Eq. d36| ) is always smaller than the tem- 
perature). Another important energy scale is the finite 
magnetic field B that is used in the simulations (i.e. how 
close the system is to half filling). We typically used a 
value of B = 0.003 which makes the magnetic length 
scale in Eq. ([To]) always negligible compared to the finite 
temperature correlation length. 



V. CONCLUSION 

We have considered a number of impurity models 
and were able to extract detailed information about 
the backscattering amplitude, the backscattering phase- 
shift, and the impurity screening effects by examining 
the Friedel oscillations. The results for the various im- 
purities have direct and indirect implications for a large 
number of theoretical models and experimental systems 
as we will summarize below. 



A. Kondo-type impurities 

Kondo impurity problems are maybe the most famous 
examples of impurity renoxmalization effects ever since 
the classic work by Wilson.EH Many of the impurity mod- 
els we have considered here are analogous to Kondo im- 
purity problems in terms of the field theory language. 
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In particular, the field theory description of a Heisenberg 
chain is the same as that of the spin-channel for a spin-full 
electron field (while the charge excitations are neglected). 
Moreover, it is known that coupling the open end of a 
Heisenberg chain to an impurity spin produces the_same 
impurity operators as in the real Kondo problem.DEd The 
number of channels in the equivalent Kondo problems is 
given by the open ends that the impurity spin is con- 
nected to (e.g. the two link impurity in Sec. IIIB is anal- 
ogous to the two channel S=l/2 Kondo problem). It 
is important to realize that the Heisenberg spins in the 
chains that we consider here have different expressions 
in terms of the boson fields than the real electron spins 
in the full three dimensional Kondo problems. Nonethe- 
less, we can still use our models to gain some insight into 
the central aspects of renormalization, scaling, cross-over 
temperature, and screening clouds. 

We have shown that the Kondo-type impurities indeed 
show the expected renormalization to a screened impu- 
rity spin. In particular, we have found a diverging screen- 
ing cloud (and vanishing backscattering) for the over- 
screened case in Sec. [II B|, while the exactly screened 
III El are characterized by a finite 



m S ec 
cases in Sec. [Ill C| and 



screening cloud and a phase shift in the backscattering 
as T -> 0. 

To analyze the renormalization process more quantita- 
tively it is important to introduce the concept of scaling. 
It can be expected that the impurity introduces a new 
energy scale that depends on the initial bare coupling 
constants. Commonly this energy scale is referred to as 
the cross-over temperature Tk- By making use of scale 
invariance it is then possible to describe the renormal- 
ization process universally in terms of the single param- 
eter T/Tk- In particular, impurity properties like the 
impurity susceptibility are described by a universal scal- 
ing function %i m p = f(T/T K )/T, which is valid for all 
T and Tk below the cut-off. This behavior was demoru 
strated explicitly before for the two weak link problems 
and works for all Kondo-type impurities in this paper 
(not shown). In fact it is possible to extract the cross- 
over temperature Tk up to an arbitrary overall scale ex- 
plicitly hjt-collapsing the data according to the scaling 
analysis old We have determined Tk this way as a func- 
tion of coupling J' in each case as shown in Fig. ^l] (up 
to an arbitrary overall scale). The Kondo temperature 
shows the same exponential dependence for small J' 



T K oc exp(-0.85J/J') 



(37) 



as shown in Fig. |2l] (coming from the same marginally 
relevant operator at the unstable fixed point in all cases). 
The underscreened case of a spin-1 coupled to the end of 
one chain has also been included in Fig. |2l| for complete- 
ness. 
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FIG. 21. Crossover temperature Tk of four different 
Kondo-type impurities. Tk has been multiplied by arbitrary 
constants in order to compare the four cases. 

More interesting in the context of the density oscilla- 
tions is maybe the scaling of the screening cloud. As the 
screening cloud we define that part of the alternating den- 
sity that is induced by the magnetic impurity, labeled by 
ci in Eqs. (§|), @, and (§f). In Ref. | it was postulated 
that the screening cloud in the real Kondo effect should 
be a function of the scaling variables xT and T /Tk- In 
our cases we can make a similar argument except that we 
need to include an overall factor T 9 ^ 1 to account for the 
dimensionality of the correlation functions. We therefore 
obtain the following scaling law 



X 



screening 



T^fixT, T/T K ). 



(38) 



Indeed we find that the shape of the screening cloud is 
not affected by Tk and can always be expressed as a 
function of the scaling variable xT. The coefficient cj 
must therefore be a function of T /Tk multiplied by ap- 
propriate powers of T . As an example we can take the 
two link problem at g = 1/2 with the screening cloud 
given in Eq. (p3|), where the coefficient can be written 
as cj = f(T K /T)/VT with some function /. In Fig. || 
we replot the coefficient cj analogous to Fig. || but with 
the argument replaced by Tk/T instead of J'. The in- 
set shows that the data indeed collapses if multiplied by 
VT as implied by Eq. (|3|) . The solid line in Fig. || there- 
fore is proportional to 1/ \JTk and diverges exponentially 
with J' according to Eq. (|37]) . Similar arguments can be 
made for the coefficients cj in the screening clouds of the 
exactly screened cases in Eqs. ( p5| ) and (p8|), except that 
cj = J(Tk/T)/T and the solid line is proportional to 
1/Tk in that case. 
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FIG. 22. The coefficient a for the two link impurity in 
Eq. (^3j) as a function of Tk /T for different temperatures 
T/J = 0.2, 0.1, 0.04, 0.025, 0.0167, 0.01 from below. 



B. Doping in spin chains 

Our results also have immediate experimental conse- 
quences for impurities in spin-chain compounds such as 
KCUF3 or Sr2CuC>3. The spin density oscillations are di- 
rectly linked to the local Knight shifts (susceptibilities) 
close to the corresponding impurities, which can be mea- 
sured by standard NMR techniques or muon spin reso- 
nance. NMR experiments have already successfully de- 
tected the sharp feature corresponding to the maximum 
in Fig. ^ from open boundaries due to non-magnetic de- 
fects that were naturally present in the crystals We now 
propose to use intentional doping with magnetic or non- 
magnetic impurities to see the predicted renormalization 
effects. Impurities of one or two modified links in the 
chain can possibly be created by doping the surrounding 
non-magnetic atoms in the crystal at link or site p arity 
symmetric locations. The spin-1 impurities in Sec. HIE 



could be produced in a more straightforward way by sub- 
stituting Cu ions by N i ions in the corresponding com- 
pounds. In Sec. 



HIE 



we discussed explicitly how the 
renormalization effects for spin-1 impurities would show 
up in an actual experiment. Similar arguments can also 
be made for the two linkEl or one link impurities by simply 
using the analytic form of the corresponding alternating 
spin densities with the coefficients cr and c/ that we have 
calculated. 

In general we find a strong enhancement of the antifer- 
romagnetic order near impurities. This enhancement can 
also be observed in higher dimensionaj and may have im- 
portant consequences for impurity-impurity interactions. 
In one dimension this effect is strongest, but the complex 
functional dependence we found here is often beyond the 
intuitive explanation in terms of valence bond states.0 



C. Impurities in Mesoscopic systems 

Finally, our analysis also allows us to draw impor- 
tant conclusions for transport measurements in one di- 
mensional mesoscopic structures. This is probably the 
first time that the conductivity could be explicitly ex- 
tracted from numerical data for Luttinger Liquid type 
models. Not surprisingly, we found that a generic im- 
purity indeed renormalizes to complete backscattering as 
the temperature is lowered, and we also could explicitly 
observe the "healing effect" in the symmetric pjesonant 
tunneling case as predicted by Kane and Fisher .Eltiti Our 
numerical results not only confirm the asymptotic power- 
laws, but also give a quantitative estimate of the conduc- 
tivity for all temperatures and impurity strengths. For 
a generic impurity with little or intermediate backscat- 
tering we hnd that the asymptotic scaling region turns 
out to be extremely narrow. For impurities with strong 
backscattering we find that the conductivity is enhanced 
by interactions at higher temperatures. 

One obvious question is how those results can be gener- 
alized to spinful electron systems and carbon nanotubes. 



A number of works ha 
purities in spinful wire: 



cessed the question of im- 
3 and found a richer struc- 



ture since renormalization takes place in both the spin 
and the charge channels. However, if realistic SU(2) in- 
variant interactions are assumed the generic behavior is 
very similar to the spinless case, so that we expect that 
our results for the reflection coefficient carry over in a 
straight forward fashion. The shape and amplitude of 
the density oscillations, however, will in general be very 
different for spinful electron systems. For carbon nan- 
otubes it has been shown that the Friedel oscillations 
impose a characteristic pattern that can be observed with 
scanning tunneling microscopyO For spinful wires it is 
expected that the Friedel oscillations from an open end 
can reveaL-the nature of the spin-charge separation in 
real spaceJ23 Although our results do not allow for quan- 
titative predictions of the density oscillations in spinful 
systems, we generally expect that strong, long-range den- 
sity oscillations should be present from backscattering in 
one-dimension. One experimental consequence of those 
oscillations is that the measurement through a lead close 
to an impurity is very sensitive to the exact location. 
Previous studies have shown that even the distance be- 
tween two leads can play a crucial roleO The current 
may be strongly enhanced or depleted, depending on if 
the distance to the impurity is a multiple of 2k px or not. 
Especially interesting arc therefore experiments with an 
adjustable lead such as a tunneling tip. The direct ob- 
servation of those oscillations could give detailed infor- 
mation about both the nature of the impurity and also 
about the interactions in the system. 
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